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Abstract 

If spacetime possesses extra dimensions of size and curvature radii much larger than 
the Planck or string scales, the dynamics of these extra dimensions should be governed 
by classical general relativity. We argue that in general relativity, it is nontrivial to 
obtain solutions where the extra dimensions are static and are dynamically stable to 
small perturbations. We also illustrate that intuition on equilibrium and stability built 
up from non-gravitational physics can be highly misleading. For all static, homogeneous 
solutions satisfying the null energy condition, we show that the Ricci curvature of 
space must be nonnegative in all directions. Much of our analysis focuses on a class 
of spacetime models where space consists of a product of homogeneous and isotropic 
geometries. A dimensional reduction of these models is performed, and their stability 
to perturbations that preserve the spatial symmetries is analyzed. We conclude that 
the only physically realistic examples of classically stabilized large extra dimensions 
are those in which the extra-dimensional manifold is positively curved. 
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1 Introduction 



The idea that spacetime may have more than 4 dimensions has had a long and distinguished 
history. In the context of string theory or other quantum theories of gravity, the extra 
dimensions typically would be expected to be of a size set by the string or Planck length. In 
such cases, one would not expect a classical description of the extra dimensions to be valid. 
However, in recent years, the possibility has been raised that "large" extra dimensions may 
be present, and a number of such explicit models have been proposed see p[ for a review, 
and for recent experimental limits). It seems reasonable to expect that classical general 
relativity should be adequate to analyze the equilibrium and stability of these models. 

In general relativity, it is nontrivial to obtain models with static extra dimensions, and it is 
even more nontrivial to stabilize these extra dimensions with respect to small perturbations. 
Furthermore, intuition developed from non-gravitational physics with regard to equilibrium 
and stability can be highly misleading. Indeed, the difficulties in obtaining stable compact 
dimensions and the counter-intuitive behavior of many models are well illustrated by a simple 
example. 

Consider a static, flat, ra-torus universe. (In this example, the n-torus will be taken to 

be all of space, but a similar example could be given where the n-torus comprises only the 

extra dimensions.) This spacetime is, of course, a static, equilibrium solution of the vacuum 

Einstein equation. Is this solution stable? The following argument suggests that it should 

not be: There is no scale that sets the values of the "moduli" of the torus. In particular, 

there is an equilibrium solution for all possible values of the radii of the torus. Therefore, 

at best, one would expect the torus to be neutrally (un)stable: If one of the dimensions of 

the torus is perturbed toward collapse, there should be nothing to prevent that dimension 

from collapsing to zero size in a finite time. It is not difficult to prove that this expectation 

is correct: There exist exact solutions of the vacuum Einstein's equation corresponding to a 

linear decrease (or growth) of the size of one of the dimensions of the torus with time,|^ so 

an arbitrarily small perturbation will cause the torus to collapse within a finite time. 

^ These exact solutions with hnear expansion or contraction of one of the radii of the torus are actuaUy 
flat. They can be obtained by starting with a wedge of (1 + 1) dimensional Minkowski spacetime where 
the action of Lorentz boosts is spacelike (the "Milne universe" ) and making a periodic identification under 
Lorentz boosts. 
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Can the above n-torus universe be stabilized via the addition of suitable matter fields? 
One might expect that it would be straightforward to do so by simply adding some matter — 
such as iron struts — whose total energy was minimized at a particular chosen set of values 
of the radii of the torus. This matter would then set a suitable equilibrium scale for the 
torus, and one might expect that the ra-torus universe at the energy minimizing radii would 
be a stable, equilibrium configuration. (In the somewhat different context of a single extra 
dimension with orbifold boundary conditions, this is essentially the idea of the Goldberger- 
Wise stabilization mechanism 0].) However, it is easy to see that this mechanism does not 
work in this case. In fact, we will see from eq. (|^) below that if one adds matter satisfying 
the strong energy condition with strict inequality holding at at least one point, then no static 
solutions whatsoever can exist, so all solutions must be expanding or contracting in at least 
one spatial dimension. Furthermore, if we consider a contracting solution with the same rate 
of contraction in all directions (thereby corresponding to a fiat FRW solution with periodic 
identifications), then it also can be seen that the presence of a positive pressure enhances 
the rate of collapse in the sense that, with the same initial conditions, a universe filled with 
matter possessing positive pressure will collapse to zero size faster than a universe filled with 
dust. Thus, rather than stabilizing the torus, adding "iron struts" or other matter with 
positive pressure actually enhances the tendency to collapse. 

The difficulties in achieving stable equilibrium of compact dimensions are not special 
to the above n-torus example. Indeed, a simple constraint on the existence of stationary 
solutions can be obtained as follows. Consider a stationary spacetime, with timelike Killing 
field t", which is spatially compact or compactifiable (i.e., possessing appropriate spatial 
symmetries so that the noncompact directions could be compactified via identifications). 
This encompasses many models considered in the literature, including all of the ones we 
shall consider below in this paper. Then satisfies (see, e.g.. Appendix C of 0) 



Since V atb = V[at6], this equation may be rewritten in differential forms notation as 



where {R ■ t)a = Rabt'^- Integrating this equation over a compact (or compactified) spacelike 




d * dt 



2*{R-t) 



(2) 
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slice S, we obtain 

J^*(R.t) = (3) 

i.e., 

/ RabU^t^dJ: = (4) 

where n°- denotes the unit normal to S, and denotes the volume element on S induced 
from the spatial metric on E. In particular, if the spacetime is static (rather than merely 
stationary) and we choose S to be orthogonal to t", then we obtain 

i^afen'^nVdS = 0, (5) 

where V = {-tHaf^'^. 

Eq. and/or eq. (P) provide a significant integral constraint on the existence of station- 
ary solutions with compact or compactifiable spacelike slices. The strong energy condition 
is the requirement that [Tab — \Tgab)v'^v^ > for all timelike v"". It is violated by a positive 
cosmological constant and by matter fields whose stress-energy is dominated by (positive) 
potential energy, but is satisfied by most other forms of matter usually considered. In the 
presence of Einstein's equation, the strong energy condition is equivalent to the condition 
that Rabv"'v^ > for all timelike f". Eqs. (^) and (^ indicate that it will be difficult to 
produce stationary solutions to Einstein's equation with matter without violating the strong 
energy condition. In particular, as already mentioned above, in the static case, if there is 
a point where Rabn°"n'' > 0, then there also must be a point where RabU^'n'^ < 0, thereby 
violating the strong energy condition. 

Similar conclusions concerning the difficulties with the existence of static solutions satis- 
fying the strong energy condition can be drawn by appealing to the singularity theorems. In 
particular. Theorem 9.5.4 of asserts that if a spacetime possesses a compact spatial slice, 
if matter satisfies the strong energy condition, and if, in addition, there are no closed timelike 
curves and the timelike and null generic conditions hold, then the spacetime must be singular 
in the sense of timelike or null geodesic incompleteness. Consequently, all stationary globally 
hyperbolic spacetimes possessing a compact or compactifiable Cauchy surface — such as the 
static ra-torus discussed above — must violate the generic condition. By the above theorem, 
any perturbation of such a spacetime that satisfies the generic condition must produce a 
singularity. 
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The above arguments against the existence of stable, static solutions can be evaded 
by allowing for the presence of a positive cosmological constant or other forms of matter 
violating the strong energy condition. Even so, it is still far from obvious how to make 
stable, static models. Indeed, the Einstein static universe provides a good illustration of 
the kind of difficulties one can encounter: Although a static solution can be obtained by 
"balancing" a positive cosmological with an ordinary form of matter, such an equilibrium is 
unstable. 

The main purpose of this paper is to investigate the constraints placed by classical general 
relativity on the existence and stability of static solutions to Einstein's equation in the case 
where "space" is taken to be a product of an arbitrary number of homogeneous and isotropic 
geometries of arbitrary curvature — in other words, conventional Kaluza-Klein theory 0. In 
sections 2 and 3, we restrict attention to the case of a product of two such geometries — one 
of which can be taken to be ordinary, macroscopic space, and the other of which can be taken 
to be the "large" (compared with the Planck scale) extra dimensions. In section 2, we show 
that if matter satisfies the null energy condition {Tabl'^l^ > for all null Z"), then existence 
of a static solution requires either that the compact space be positively curved, or that 
the spacetime be locally Minkowski and the stress-energy tensor vanish identically. Thus, 
provided only that the null energy condition holds, the extra dimensions in the static models 
with the symmetries we consider cannot have negative curvature, and they cannot even be 
fiat unless the universe is completely empty. More generally, we show that for any static, 
homogeneous solution with matter satisfying the null energy condition, the Ricci curvature 
of space must always be nonnegative. 

In section 3, we consider a particular model that has previously been considered in the 
literature |^, ^ |^, [10], |Tl]], in which spherical extra dimensions are stabilized via a balance 
between vacuum energy and gauge fields (and the curvature of the internal space). We 
perform a dimensional reduction of it, thereby reducing it to an ordinary Robertson- Walker 
model with a self-interacting scalar field. The static solution is unstable in the sense that 
an arbitrarily small perturbation of it will result in a ( "big bang" ) singularity. However, the 
perturbed model appears to be cosmologically acceptable. 

Section 4 is devoted to exploring the relationship between energy minimization versus 
equilibrium and stability in general relativity. In non-gravitational physics, a configuration 
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which minimizes the total matter energy corresponds to a stable equilibrium configuration, 
and stability arguments are typically made on the basis of energy minimization. How- 
ever, such arguments are, at best, highly suspect in general relativity. In asymptotically 
fiat spacetimes — where total energy is well defined — configurations of minimum total en- 



ergy should correspond to stable equilibrium configurations |jT2[ , but the total energy of the 
spacetime is not equal to the integrated energy of matter in the spacetime, so minimization 
of the total matter energy may have little to do with the conditions for stable equilibrium. 
The situation is considerably worse in compact or compactifiable spaces, where no nontrivial 
notion of total energy exists. To investigate the relationship between minimization of total 
matter energy and stable equilibrium, we consider a Robertson- Walker spacetime filled with 
a fiuid, and ask the following two questions: (i) What conditions on the pressure of the 
fiuid must be satisfied in order that the total energy of the fiuid be minimized? (ii) What 
conditions on the pressure of the fiuid must be satisfied in order to have a stable, static 
solution to Einstein's equation? We shall see in section 4 that, in general, the conditions 
for equilibrium in general relativity bear no relationship with the conditions for energy ex- 
tremization. Interestingly, however, in 2 + 1 dimensions, it turns out that the conditions for 
equilibrium in general relativity actually coincide with the conditions for energy extremiza- 
tion. Remarkably, in this case we find that the conditions for stable equilibrium in general 
relativity are precisely that the matter energy be a local maximum — exactly the opposite of 
what would be expected in non-gravitational physics. 

Finally, in section 5, we generalize the results of sections 2 and 3 to the case where 
the spatial geometry is a product of an arbitrary number of homogeneous and isotropic 
geometries. We show that, again, if matter satisfies the null energy condition, the existence 
of a static solution requires nonnegative curvature of each of the geometries; it is not possible 
to compensate for some negative curvature by the addition of positive curvature in a different 
factor. Following |T^, |I^ we also perform a dimensional reduction in this case, re-expressing 
the model as an ordinary Robertson- Walker cosmology with a collection of self-interacting 
scalar fields. Considering a specific example with four extra dimensions in the form of 
a product of two two-spheres [^, we investigate the cosmological acceptability of such 
compactifications. 

Our analysis leaves numerous questions unanswered. We do not, for example, consider 
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off-diagonal perturbations to our product spacetimes or perturbations that do not respect the 
spatial symmetries of the extra dimensions. Also, models where the standard-model fields 
are confined to a "brane" with a thickness much smaller than the size of the extra dimensions 
are incompatible with our assumption of perfect homogeneity in the extra dimensions. (In 
particular, Randall- Sundrum spacetimes with a single finite extra dimension rely on the 
tension of the branes on the boundary, and are therefore not addressed by our analysis.) 
While a natural first step is to assume that the energy localized on the brane is smaller than 
that in the "bulk" and can hopefully be neglected, it is important to subsequently check 
that the introduction of a nonzero brane tension does not destabilize the solution. (Some 
progress in including the effects of nonzero brane tension has been made by considering 
metrics generated by global topological defects |I5[].) 



2 Static, homogeneous solutions in general relativity 

In this section we examine static homogeneous spacetimes, without regard to questions of 
stability (considered in later sections). We consider spacetimes of the form H. x Aii x M.2-, 
where R is the timelike direction and M.i and M.2 are maximally symmetric manifolds of 
dimensionality N and n, respectively. We can think of A^i as representing the "large" spatial 
dimensions and M.2 as representing the "extra" ones, although they will be treated on an 
identical footing. The metric can be written 

ds^ = —dt^ + gijdx^dx^ + jijdy^dy^ , (6) 

where gij{x^) and '^ij{y^) are metrics on manifolds of constant curvature (either positive, 
negative, or zero). Note in particular that gij{x^) and 'Jijiy^) are independent of time. 
Indices /, J run from 1 to A^, while i,j run from 1 to n. We define curvature parameters by 

^ Rjgij] ^ ^[7»i] .7^ 

N{N-1)' n{n-iy 

where R is the Ricci scalar constructed from the appropriate spatial metric. 

For a stress-energy source we take a perfect fluid defined by an energy density p and 
pressure in the large and small dimensions, P and p: 

Too = P (8) 
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Tij = Pgij (9) 

Tij = p-fij. (10) 

There are three independent nontrivial components of Einstein's equation, corresponding to 
the 00, IJ and ij components: 

^N{N -1)K +^n{n-l)k = 8nGp (11) 

^{N -1){N -2)K +^n{n-l)k = -StiGP (12) 

^N{N -1)K + ^{n-l){n-2)k = -SnGp, (13) 



where G is Newton's constant (in A^ + n+ 1 dimensions). In ([TT|-[T3|) we have put overbars on 
the energy density and pressures to emphasize that these relations hold at a static solution. 

We next impose the null energy condition^] (NEC), that Tahiti'' > for all null vectors l". 
In the context of our model, the NEC is equivalent to the requirements 

p + P > 0, (14) 
p + p > 0. (15) 

Taking appropriate linear combinations of (|ll]-0) yields 



8nG{p + P) = {N-1)K, (16) 
8tiG{p + p) = {n-l)k. (17) 

From comparing (|TB|-|I^ with (|T^-[T^), we see immediately that the curvatures K, k must 
be either zero or positive. (The curvature cannot be made negative by choosing a one- 
dimensional manifold, as a one-dimensional space must have zero curvature.) Thus, in a 
static spacetime obeying the NEC, with spatial sections described by a product of two 
symmetric spaces, neither spatial factor can be negatively curved. In particular, the NEC 
must be violated in order to consider stabilization of compact hyperbolic extra dimensions. 



which have recently been investigated [16 



^The NEC is in fact weaker than the "weak" energy condition, that TabV°'v^ > for all timelike vectors 
v"" . By continuity the weak energy condition implies the NEC, but not vice-versa. In particular, the NEC 
permits a negative energy density if there is an equal and opposite pressure, and thus allows for a cosmological 
constant of either sign, while the weak energy condition prohibits any negative energy density, and thus allows 
only a nonnegative cosmological constant. 
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Several possibilities remain open: both spaces may be fiat, one may be fiat and the other 
positively curved, or both may be positively curved. If both spaces are fiat {K = = 0), 



we see from (|TT]-|1J) that all of the components of the stress-energy tensor must vanish, 
p = P = p = 0. This solution is (locally) empty Minkowski space. 

The nonempty static solutions obeying the NEC are thus those in which at least one 
space is positively curved, and neither is negatively curved. This implies that toroidal com- 
pactification [k = 0) is only possible if the large dimensions are positively curved {K > 0). 
In fact, however, this case is not physically realistic, for two reasons. First, in the real uni- 
verse the three observed spatial dimensions are fiat to an excellent approximation. (Evidence 
for spatial flatness has been provided by measurements of the anisotropy spectrum of the 
cosmic microwave background [|17|; for our present purposes, however, all we need to know 
is that radius of curvature of the large dimensions is significantly larger than any relevant 
microphysical scale.) Second, when (in the next section) we consider stability of a model 
with both flat and positively-curved dimensions, we will find that small perturbations induce 
oscillations in the positively curved space and a big-bang singularity in the flat dimensions. 
Hence, a toroidal set of extra dimensions would be unstable to collapsing to zero volume 
(or expanding to infinite volume). We therefore conclude that the only physically relevant 
examples of large extra dimensions obeying the NEC are those with positive curvature. 

The above result on the non-negativity of spatial curvature for static solutions can be 
significantly generalized as follows. Consider any static, spatially homogeneous spacetime, 
with static Killing field t". Now, any static Killing field satisfies 

Vah = -2V-%V,]V (18) 

where V = (—t°'taY^'^. However, for a spatially homogeneous spacetime, we have = 0, 
so t°- is covariantly constant. It then follows from eq. (^ that Rabt'' = 0, where Rat denotes 
the Ricci curvature tensor of the spacetime metric. Let be any unit "spatial" vector (i.e., 
= 0), and let l" = V^^t" + s". Then Rabl^'l'' = RabS^sK However, since I'' is null, by 
Einstein's equation, we have 

Rabtl' = SnGiTab - \gabT)rf = SnGTabtl' (19) 

Consequently, the NEC together with Einstein's equation implies RabS°'s^ > for all spatial 
vectors s"", with equality holding if and only if Tabl°'f = 0. Furthermore, since is covariantly 
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constant, we have Rabs'^s'^ = Rabs"'s^, where Rab denotes the Ricci curvature of space. Thus, 
for any static, spatially homogeneous spacetime satisfying the null energy condition, the 
Ricci curvature of space must always be nonnegative. 

3 An example: stabilizing 5^ with gauge fields 

Having established the criteria for static solutions with a single extra-dimensional space, we 
now will examine in detail a specific example to understand the sense in which it can be 
"stable" . Our example will be a spacetime with four macroscopic spacetime dimensions and 
two positively curved extra dimensions, stabilized by balancing a cosmological constant and 
a gauge field wrapped around the extra dimensions. (This mechanism was recently exam- 
ined by Sundrum [|11[]; a nonabelian version was originally investigated by Horvath, Palla, 
Cremmer and Scherk |^, which was generalized to larger numbers of extra dimensions by 
Freund and Rubin Q], and stability was discussed by Randjbar-Daemi, Salam and Strathdee 
1^. A thorough discussion, including multiple compact spaces, is given by Guenther and 
Zhuk [|TU|.) We consider the conventional dimensional reduction of this theory, rewriting the 
six-dimensional Lagrangian as a four- dimensional Lagrangian plus a scalar field (the radion) 
representing the size of the extra dimensions. The effective potential for the radion will be 
found to have two minima: one at a finite radius, and one at infinite radius. This kind of 
analysis is by no means original (see for example [0 and references therein), but we will go 
through it carefully for pedagogical purposes. 

Let Gab be the metric for the full spacetime with coordinates X", and denote by Gn+u+i 
Newton's constant in the full spacetime. We consider N large spatial dimensions and n extra 
dimensions, so that indices a, b run from to + ra; we will specialize to the case = 3, 
n = 2 after performing the dimensional reduction. We consider metrics of the form 

ds"^ = GabdX^dX^ = gAB{x)dx^dx^ + b'^{x)^ij{y)dy'dy\ (20) 

where the are coordinates in the (A^ + l)-dimensional spacetime and the are coordinates 
on the extra-dimensional manifold, again taken to be a maximally symmetric space with 
metric 7jj. However, for the dimensional reduction analysis given below, it will not be 
assumed that the metric qab of the large dimensions has any symmetries. The action is the 
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(A^ + n + l)-dimensional Hilbert action, plus a matter term we leave unspecified for now: 



S = / f/^+"+iXv^ — — R[Gab] + Cm , (21) 

where \/—G is minus the square root of the determinant of Gab-, R[Gab\ is the Ricci scalar of 
Gab) and Cm is the matter Lagrangian density with the metric determinant factored out. 

The first step is to dimensionally reduce the action (PT]): we express everything in terms 
of qabi liji and 6(x), and integrate over the extra dimensions 0. From the metric ansatz 
(pop we have 

= fe"v^V7, (22) 
R[Gab] = R[gAB\+h-^R[^,,]-2nh-^g^''VAVch-n{n~l)h-^g''''{VAh){Vch), (23) 

where Va is the covariant derivative associated with the (A^ + l)-dimensional metric qab- 
We denote by V the volume of the extra dimensions when 6 = 1; it is given by 



V = / rf"2/V7. (24) 



The (A^+ l)-dimensional Newton's constant Gn+i is related to its higher-dimensional version 
by 



IGvrGTV+i IGvrG 



V 

(25) 



N+n+l 



We are thus left with 



J 167rG;vj_i L 



IGTrGAT+i 

-n{n - 1)6"- V^(Va6)(Vb6) + n{n - l)fc6"-2] + V6"£a/|, (26) 

where the curvature parameter k of 'jij is given by (|^). 

■^As is well known, if one substitutes a metric (or other field) ansatz into an action, there is a danger 
that variation of the new action will not reproduce all of the field equations, since the metric (or other field) 
variations are restricted by the ansatz. If the ansatz arises from symmetry assumptions, then, typically, the 
missing field equations will be automatically satisfied by virtue of the assumed symmetry. It is not difficult 
to verify that this is the case here, so that the reduced action that we will obtain below does indeed yield 
all of Einstein's equations. On the other hand, if the ansatz arises, e.g., from a gauge choice, then there is 
no reason to expect that the new action will reproduce all of the field equations. An example of this can be 
found in the last reference of 0|, where the gauge choice goo = 1 is imposed in the action, and the action 
consequently fails to yield the 00-component of Einstein's equation. 
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Through a change of variables and a conformal transformation, 

b{x) = e^(^), (27) 
9AB = e-^-^l^'^-^AB. (28) 

we can turn this reduced action into that of a scalar field coupled to gravity in the so-called 
"Einstein frame" (where the gravitational Lagrangian is simply the Ricci scalar, with no 
multiphcative scalar fields) . Some algebra yields 

+n(n - \)ke-'^^^^^\ + Ve-'^^^^l^, (29) 

where we have dropped terms that are total derivatives. 

To turn /? into a canonically normalized scalar field, we make one final change of variables, 

to 

^=^ "'^/_7^> Mpft (30) 

where the reduced Planck mass is Mp^ — ^J^'kG jv+i (in units where h — V). We are then 
left with 



S = yd^+ixy^|^M|i?[5^s]-i^^^(VA</>)(VB0) 



2 



Specializing to the case n = 2 and = 3, we see that is a canonically normalized scalar 
with a potential 

V{<\>) = -kMj.e-^^/^'^ - Ve-'^/^^£M(0). (32) 

The stabihty properties and cosmological evolution of this model will evidently depend on 
the choice of matter Lagrangian jCm{(P)- 

As an attempt at a realistic scenario, we take Cm to consist of a cosmological constant 
plus an abelian gauge field wrapped around the extra dimensions. In the six-dimensional 

picture, we have an antisymmetric field strength Fat which we take to be proportional to the 
volume form of the extra-dimensional metric ; in components we have 

F45 = -i^54 = (33) 

12 



where Fq is a constant and all other components vanish. In particular, the components of 
Fab are independent of b. (It is straightforward to check that such a field satisfies Maxwell's 
equations, either directly or by using conservation of flux.) We have a six-dimensional matter 
Lagrangian 

£m = —^6 — -^G'^'^G^'^FabFcd, (34) 

where Ag is the six-dimensional vacuum energy density. Upon dimensional reduction this 
yields 

V£m = -A4-/oV2^/*^-, (35) 

where A4 = VAg is the "bare" four-dimensional vacuum energy density and /g = \VFq is a 
(nonnegative) constant. 

The effective potential for in this model is therefore 

r (0) = - A;M|e-2<^/*^^ + X^e-'^""^ + fle-^^"'^ . (36) 

We would like to choose the parameters A4 and /o such that there is a local minimum of 
the potential for which = (so that the effective four- dimensional cosmological constant 
vanishes and we can have a fiat-spacetime solution). It is clear from ( PBD that choosing k < 
cannot lead to an isolated minimum at finite (p, in accordance with the results of the previous 
section. Instead we choose k > and 

A4 = = \kMl. (37) 

There is then a local minimum at = 0, and the potential takes the form shown in Figure |l|; 
the mass of the radion is 

ml = ^ , _ = k. (38) 



<j>=0 



Note that there is also a minimum as (j) goes to infinity, corresponding to blowing up the 
size of the extra dimensions. This raises the possibility of a Brustein-Steinhardt problem [|19|, 
if cosmological evolution tends to overshoot the desired minimum at = 0. However, this 
does not seem like a serious obstacle to the most popular scenarios, which feature two 
extra dimensions, a six- dimensional Planck scale of order 1 TeV, and a volume V of order 
(1 mm)^ ~ (10^^ eV)"^, so that the four-dimensional Planck scale is Mp ~ 10^^ GeV [0. 
From (PSP the radion is a light field, with mass set by the length scale of the extra dimensions. 
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Figure 1: The effective potential ( ^6]) of the radion in the model described in the text. 

Meanwhile, the height of the barrier in Figure [l| is of order 1^(0 = Mp) ~ kMp ~ 1 TeV^. 
Thus, in order for the field to roll coherently over the barrier, it would have to have an 
energy density of order the (six-dimensional) scale of quantum gravity in this model, which 
is in a regime in which our classical analysis has no right to be applicable. However, another 
potential difficulty with this model arises from the possibility that perturbations in the early 
universe could lead to localized regions in which tunnels over the barrier, either due to 
thermal or quantum fluctuations, giving rise to "bubbles" in which the extra dimensions 
expand to infinity. We do not analyze this issue here. 

Thus, with this model, the effect of the extra dimensions is to introduce an additional 
massive scalar field (the radion) into the four- dimensional theory. It is therefore straightfor- 
ward to answer questions concerning stability; if we consider a solution in which the radion 
is slightly perturbed from its equilibrium value, it will act as an energy source in the large 
dimensions, which will evolve as a Friedmann-Robertson- Walker spacetime. This spacetime 
will have a Big-Bang singularity, but is nevertheless cosmologically acceptable. The radion 
acts as a massive scalar which contributes to the total energy density of the universe; so 
long as its contribution is not too large (or it can decay away), it will not be cosmologically 
harmful. 
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4 Stabilization vs. minimization of energy 



When one is attempting to stabilize the extra dimensions by adding some form of matter, 
it might appear natural to expect that the extra dimensions will be stable when the total 
matter energy is minimized (the total matter energy being the integral of the energy density 
of matter over the volume). Indeed, this conclusion is rigorously correct for systems normally 
considered in non-gravitational physics: the total energy of matter normally is positive 
definite (thereby providing an appropriate norm) and conserved, and only positive energy 
can be radiated to infinity. Consequently a small perturbation about the state of minimum 
energy cannot drive the system far from this state. However, this argument fails in general 
relativity, since there is no global conservation law for the total energy of matter^. The 
main purpose of this section is to consider a simple cosmological model illustrating how one 
can be led seriously astray by non-gravitational intuition when determining the stability of 
spacetimes in general relativity. 

Consider an (n + l)-dimensional spacetime of constant spatial curvature k, so that the 
metric takes the conventional Robertson- Walker form 

ds'^ = -dt^ + a{tfja(3dx"dx'^. (39) 

This space is to be filled with a uniform, perfect fluid of energy density p, whose pressure, 
p, is a function of p. By conservation of stress-energy, we have 

^dp 

P=-P--^. (40) 
n da 

This equation can be integrated to express p (and, hence, p) in terms of a, and, in the 
following, it will be convenient to view p and p as functions of a. The quantity U = pa^ 
represents the "total energy" (per comoving volume element) of the fluid. Minimizing U 
with respect to a, we obtain 

^ = -pna"-\ (41) 
da 

Thus for an extremum of U we must have p = 0. Furthermore, for this extremum to be a 
minimum, we require (PU /da^ > 0, giving 

d'^U „ -idp , , , 

= -na^'-^-f- > 0, (42) 
p=o da 



da^ 



^Inflation provides a dramatic illustration of this fact. 
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which of course is just the condition that dp/ da be negative. Thus, analysis of the total 
matter energy in the manner normally done in non-gravitational physics suggests that space 
can be in equilibrium at some size only when the pressure is zero, and stable only if the 
pressure increases as the volume shrinks ("the more the system is squeezed, the harder it 
pushes back"). 

But now consider the analysis of equilibrium and stability in general relativity, as deter- 
mined by Einstein's equation. For the Robertson- Walker metric (^) we obtain the Fried- 
mann equations in (n+1) dimensions, 

-n{n-l)(-y = 87iGp--n{n~l)k, (43) 
(n-l)(-) + i(n-l)(n-2)(-)' = -87rGp-(n-l)(^-l)fc, (44) 

where k is the curvature parameter (|^ for the spatial manifold. The conditions for equilib- 
rium (d = d = and setting a = 1) are 

n — 2 

V = P- (46) 

n 

which bear little relationship to the condition p = obtained by extremization of total 
matter energy. 

To determine whether the equilibrium solution is stable we examine the linearized equa- 
tions of motion. We expand all variables to first order: a = l + 6a, p = p + 6p, and p = p + 5p. 
Then 

= 87iGSp + n{n-l)kSa, (47) 

{n-l)6a = -8nG5p + {n - l){n ~ 2)k5a 

= -8TTG^5a+{n-l){n-2)k5a. (48) 
da 

To obtain a stark contradiction to the non-gravitational analysis, set n = 2. The second of 
the equations for equilibrium (^) then reduces top = 0, which is precisely what was obtained 
by energy extremization. Therefore, we may compare the general relativistic conditions for 
stability with those that would be predicted by energy minimization. The second of the 
equations for stability ( |I5| ) gives 

5a = -8TTG^6a. (49) 
da 
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In other words, this (2 + l)-dimensional spacetime is stable only when dp/ da is positive, 
exactly the opposite of eq. (^). Thus, in this model, the conditions for stable equilibrium 
in general relativity are that the total energy of matter be a maximum. 

5 Solutions with multiple compact spaces 

In this section, we generalize the discussion of sections 2 and 3 to the case of multiple 
maximally symmetric nj-dimensional spaces. We begin by characterizing static solutions as 
in section 2. We then recast the full gravitational action as the action for a set of scalar 
fields evolving in a four-dimensional spacetime as in section 3, and discuss the possibility 
of stable solutions. Multiple compact spaces have been previously considered in [1^, [l^; to 
this discussion we add an argument that the extra dimensions should be positively curved, 
and we consider the height of the barrier separating the static solution from one in which 
the scale factor for the extra dimensions runs away to infinity. 

Our consideration of static spacetimes is a straightforward generalization of that in section 
2, so we will be concise. Consider a spacetime with spatial sections described by products 
of m maximally symmetric manifolds, such that the metric may be written 

m 

ds' = -df + J2j^}p4x^'dxf'^, (50) 

i=l 

where the 7^*^/3- 's are metrics of constant curvature. The indices run from + 1 

to rii with no = 0, i.e. the set of indices "runs over the submanifold . The curvature 
parameters ki are defined analogously to (|^). For a stress-energy source, we consider a 
perfect fluid defined by an energy density p and a pressure pi in each set of dimensions: 

Too = P (51) 
TaA=P^7Sp.■ (52) 
There are m + 1 independent components of Einstein's equation: the 00 component, 

H ^^ji^j - = ^T^Gp, (53) 
j ^ 

and the a^/^i components, 

- (n. - l)fc, + E \nAn, - l)k, = S.Gp., (54) 

17 



where once again overbars remind us that we are considering static solutions, and G is 
Newton's constant in the full spacetime. 
The null energy condition implies 



P + Pi > 0, 



(55) 



while the conditions (|53| - |54D lead to the relation 



87rG(p + Pi) = {rii - l)ki. 



(56) 



Thus, the left hand side of (0) is nonnegative, and there are no static solutions with any 



fcj < that obey the NEC, in accordance with the general argument given at the end 
of section 2. None of the spatial factors can be negatively curved; it is not possible to 
compensate for some negative curvature in one factor with additional positive curvature 
somewhere else. As before, if all of the curvature parameters vanish, we recover empty 
Minkowski space. The only static nonempty solutions will involve no negatively curved 
spaces and at least one positively curved space. 

With the criteria for static solutions in hand, we will now explore the stability of these 
solutions to small perturbations in the scale factors. Following section 3, we first consider 
the conventional dimensional reduction of this theory, rewriting the 4 + Srij dimensional 
Lagrangian as a four dimensional Lagrangian plus M scalar fields representing the size of 
each set of extra dimensions. (The sum is only over the extra dimensions, and M refers to 
the number of sets of extra dimensions.) 

Let Gab be the metric for the full spacetime with coordinates X", and denote by G4+Ar 
Newton's constant in the full spacetime where M = T,ni, the number of extra dimensions. 
The indices a, b run from to 3 + A/". We consider metrics of the form 



where A and B run over the four indices of the 3+1 dimensional spacetime and each set 
run over the appropriate submanifold. The action is as before: 



ds'^ = Gabdx'^dx^ = QABdx^dx^ + ^h1{x)'^^^]p^{y)dy'^'dy 



(57) 



i=l 



s 



I 



( 



167rG4+_^ 



1 



R[Gab] + 



] 
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From the metric ansatz we have 



M 



(59) 
(60) 



i=l 



-n,(n, - l)bfg''^idAbi)idcbi) -J2n,n,bi%'g''^dAkdcb,], 

where Va is the covariant derivative associated with gAB- The volume of the extra dimensions 
when all 6j = 1 is given by 

V = n / d'^'y (61) 



and the four- dimensional Newton's constant G4 is related to its higher- dimensional version 
by 

1 V 

167rG'4 167rG44-Ar 

We are thus left with 

Y\^bT 



M 



S = / d^x^l R[gAB] + EK(^^ - - 2mbr^g^''VAVcbi 



i=l 



- l)br'g^^{dAk){dck) -Y.n,njbr%'g^''dAhdcbj] 



+VC 



AI 



(63) 



where we have introduced the curvature parameter ki of each set of extra dimensions via 

i?[7«] = mim - l)ki. 

Through a change of variables and a conformal transformation, 

ai{x) = ln6j(x) — A (64) 
gab = e^^gab, (65) 

where A = —jUriiliabi, we can turn this reduced action into that of M scalar fields coupled 
to gravity in the Einstein frame. Ignoring a total divergence term, after some algebra the 
reduced action becomes 

1 



S 



167rG4 



(66) 
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where 

rjij = ninj/{2 + Af) - riiSij . (67) 

In order to have terms that look hke the canonical kinetic terms for a set of scalar fields, 
we must diagonalize the quadratic form rjij] this has been performed explicitly by |13 
Thus, we may rewrite the action as 

^ R[gAB] + E e-'"'n,(n, - 1)^^ (68) 



s = H-~9 



where !][ are the eigenvalues of rjij and the a[ are combinations of the a, that correspond to 
the Since each a, may be written as a linear combination of the a^, we have an action 
for M scalar fields on a background spacetime. The potential for these fields is positive for 
Cm > 0, although from this form it is not clear in general whether the static point is stable. 

As an example, consider the spacetime A4 x S"^ x S"^, where is a 3+1 dimensional 
homogeneous, isotropic model [1^. The action for this spacetime reduces to 

S = j ^g{\MlR[gAB]-\g^''dA<PdB<P-\g^''dA^dB^ 



where V is the volume of the extra dimensions as defined in eq. (|6T|) above, and we have 
introduced 

= -^(«/ + «,.)Mp (70) 
^ = {af-ah)Mp, (71) 

where / and h label the spheres. 

We choose the matter for this model to consist of a cosmological constant. As, and abelian 
gauge fields. Fab and Hab, so the matter Lagrangian is 

r \ ^/^ac/^bdrp rp ^ ryacrybdrj tt {70\ 

-L-M — — As — -tj Lt Pab-fcd—^^ Uab-tlcd, U^J 

We take the gauge fields to wrap around each S*^, so that, in components, we have 



F45 = -F54 = V7(/)Fo (73) 



Her = -ifre = V7(")^o, (74) 
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where Fq and Hq are constants, and all other components vanish. The conditions for a static 
solution are [see eqs. ( |55| , pBD ] 



SnGsp = 3K + kf + kh (75) 

87rG8(p + ^) = 2ir (76) 

SnGsip + p-f) = kf (77) 

8nGs{p + Ph) = h, (78) 

where Gg is the 8-dimensional Newton's constant, K is the curvature parameter for the 
large dimensions, and kf and k^ are the curvature parameters for jlp and \ respectively. 
Consequently, we find that static solutions exist, provided that 

A4 = ^MUkf + kh) (79) 
K = (80) 
/o = lMj,kf (81) 

hi = \Mlkh, (82) 

where 

/o = ^VFo^ (83) 

hi = \VHI. (84) 

Thus, there exists a 2-parameter family of static solutions. It is convenient to choose as 
independent parameters the higher- dimensional (reduced) Planck mass and the ratio of the 
static curvature parameters of the two-spheres: 



Mg-^ = ^MGg (85) 

(86) 



kf 



In terms of these, the equilibrium volume of the internal space is set by 
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where Mp = 10^^ GeV is the (reduced) four-dimensional Planck scale. We will assume that 
Mg = 1 TeV; this implies V ~ (30 keV)"*^. The curvature parameters for the spheres can be 
determined from (|86|) and 

kfkh 

we have for example 

kf = 4vr^v^. (89) 
The potential in this matter model is 

'l + r' 



2-kMpM. 



I ^ / -2[(2/v^)0-v]/Mp _ 2Q-{Vs<P-i')/Mp\ _ mo) 

The fields which diagonalize the mass matrix at the equilibrium solution {(p = 0, ip = 0) are 

2(r-l) 



2 - 










-2/3 







V3(r + 1)(2-^|) 
2(r - 1) 



,v^(r + l)(/i|-2/3) 



(p + ip 



(91) 



(92) 



where the /I's are normalized masses, 



1^1 



A4 

m%Ml 



3r 



(1 + r) 



; 1 



'1 + 



3r 



(93) 
(94) 



and m$ and are the actual masses of the $ and \E' fields respectively. 

It is clear that replacing r — > 1/r is equivalent to interchanging the roles of the F and 
H gauge fields; we therefore only consider r > 1, 01 kf > (so that the sphere with F- 
flux is sma//er than that with if -flux). There is also an upper limit on r, derived from the 
requirement that the curvature of the extra dimensions not exceed the higher-dimensional 



Planck scale, kf < Mg. From (p9|), setting Mg = 1 TeV then implies 

r < 10^^ 



(95) 
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There is a complementary constraint, that the larger two-sphere not be big enough that 
deviations from Newtonian gravity would have been detected in laboratory experiments; 
this corresponds to < (10~^ eV)^^ This bound is just saturated when r = 10^^, so 
these bounds are essentially equivalent. 

For finite values of r, the masses of both $ and are positive, implying stability of 
the system around the equilibrium solution. We would like to comment briefiy on the 
cosmological acceptability of such a scenario, following the discussion of the height of the 
potential barrier investigated in section 0. For a given value of r, the relevant height is that 
of the lowest barrier separating the equilibrium solution from infinity. This will correspond 
to a saddle point in V^($, \E'), and the barrier height will be the value of V{^, \E') evaluated 
at the saddle point. Numerically investigating the potential (^) with Mg = 1 TeV, we 
find that the height of the barrier varies with r, but the barrier always remains at least 
1 TeV^. Since this is the fundamental scale of quantum gravity, we again find that there is 



no Brustein-Steinhardt problem for this theory. At energy scales where classical gravity 
should be applicable, the field will never be able to classically roll over the barrier. 
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